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We analyze spatial structure of the neutron Cooper pair in superfluid low-density uniform matter 
by means of BCS calculations employing a bare force and the effective Gogny interaction. It is 
shown that the Cooper pair exhibits a strong spatial di-neutron correlation in a wide range of 
neutron density p/po ~ 10“^ — 0.5. This feature is related to the crossover behavior between the 
pairing of the weak coupling BCS type and the Bose-Einstein condensation of bound neutron pairs. 
We also show that the zero-range delta interaction can describe the spatial structure of the neutron 
Cooper pair if the density dependent interaction strength and the cut-off energy are appropriately 
chosen. Parameterizations of the density-dependent delta interaction satisfying this condition are 
discussed. 

PACS numbers: 21.60.Jz, 21.65.-|-f, 26.60.+C 


I. INTRODUCTION 


The importance of the pair correlation has been widely recognized for nucleon many-body systems in various 
circumstances, in particular in open-shell nuclei and in neutron stars. The pairing gap varies with system parameters 
such as N, Z and the rotational frequency in the case of finite nuclei, or the temperature and the density in the case 
of neutron stars (cf. Refs. dillilBl as reviews). The pairing correlation at low nucleon density is of special interest 
since the theoretical predictions for low-density uniform matter suggest that the pairing gap may take, at around 1/10 
of the normal nuclear density, a value which is considerably larger than that around the normal density iiH This 
feature is expected to have direct relevance to the properties of neutron stars, especially those associated with the inner 
crust HQ The strong pair ing at low density may be relevant also to finite nuclei, if one considers neutron-rich nuclei 
near the drip-line I^TToL lll|. This is because such nuclei often accompany unsaturated low-density distribution 
of neutrons (the neutron skin and/or the neutron halo) surrounding the nuclear surface [T^ ll.lL IT ^ . It is interesting 
to clarify how the pair correlation in these exotic nuclei is different from that in stable nuclei, reflecting the strong 
density dependence mentioned above. In this connection we would like to ask how the pair correlation at low nucleon 
density is different from that around the normal density. 

Spatial structure of the neutron Cooper pair is focused upon as a characteristic feature of the low density nucleon 
pairing. Its possible indication could be the di-neutron correlation in the two-neutron halo nu clei, e.gj . ^^Li, for which 
a spatially correlated pair formed by the halo neutrons has been predicted theoretically [p, flR IitI flM ll^ 1^ 
and debated in experimental studiesj^ iMUlil Further, a recent theoretical analysis]23 using the Hartree-Fock- 
Bogoliubov (HFB) methodp^ l28i l29| predicts also presence of similar di-neutron correlation in medium-mass 
neutron-rich nuclei where more than two weakly-bound neutrons contribute to form the neutron skin in the exterior 
of the nuclear surface. It is also possible to argue importance of the spatial correlation from a more fundamental 
viewpoint based on the nucleon interaction in the channel. The bare nucleon-nucleon interaction in this channel 
has a virtual state around zero energy characterized by the large scattering length a ~ —18 fm, which implies a very 
strong attraction between a pair of neutrons with the spin singlet configurations. A rather general argument |3(i l3l| . 
which applies to a dilute limit of any Fermion systems, indicates that the pair correlation of the Fermions interacting 
with a large scattering length differs largely from what is considered in the conventional BCS theoryassuming weak 
coupling: it is then appropriate to consider a crossover between a superfluid system of the weak-coupling BCS type 
and a Bose-Einstein condensate of spatially compact bound Fermion pairs [^|^|^|^|^. This BCS-BEC crossover 
phenomenon was recently observed in a ultra-cold atomic gas in a trap for which the interaction is controllable!^. 
In the case of the nucleon pairing, the BCS-BEC crossover has been argued mostly for the neutron-proton pairing in 
the ^SDi channel, for which the strong spatial correlation associated with the deuteron and the BEC of the deuterons 
may emerge |3T 1^ liol I4T| . Concerning the neutron pairing in the channel, which we discuss in the present 
paper, we may also expect that the strong coupling feature may lead to the spatial di-neutron correlation although 
the realization of the crossover could be marginal and depend on the density [^. 

In the present paper, we would like to clarify how the the spatial structure of the neutron Cooper pair varies 
with the density. For this purpose, we shall investigate the neutron pair correlation in symmetric nuclear matter 
and neutron matter in the low density region. Uniform matter is of course a simplification of the actual nucleon 
configurations in hnite nuclei and neutron stars. However they have a great advantage as one can solve the gap 
equation in this case without ambiguity for various interactions including the bare nucleon-nucleon forces with the 
repulsive core H, ^ |E EE EE EE EE EE EE well as effective interactions such as the Gogny 
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forcem m m, provided that the BCS approximation (equivalent to the HFB in finite nuclei) is assumed. It is 
straightforward then to determine the wave function of the neutron Cooper pair from the solution of the gap equation 
H m Ei El El El El Ei This provides us with a good reference frame to study the spatial structure of the Cooper 
pair as a function of the density while we do not intend to make precise predictions on other properties of neutron 
and symmetric nuclear matter. We shall perform an analysis employing both a bare force and the effective Gogny 
force|52j. and using a Hartree-Fock single-particle spectrum associated with the media. Our main conclusion will be 
that the spatial di-neutron correlation is strong in a wide range of the low density p/po ~ 10“^ — 0.5, independently on 
the adopted forces. We shall clarify the nature of the strong spatial di-neutron correlation in terms of the BCS-BEC 
crossover model. 

We shall also examine a possibility of phenomenological description of the spatially correlated neutron Cooper 
pair. Here we consider a contact force with ajrarametrized density dependent interaction strength, called often the 
density dependent delta interaction (DDDI) O; EE 1^ EH EH EE EH ■ The parameters of the DDDI need 

to be determined from some physical constraints. For example, the interaction str eng th has been constrained by 
conditions to reproduce the experimentally extracted pairing gap in finite nucleijE mil y, El i3, or the 
density-dependence of the neutron pairing gap in symmetric nuclear matter and the scattering length [M. . It 
should be noted here that the contact force requires a cut-off energy, which needs to be treated as an additional 
model parameter. Concerning the cut-off parameter, attentions have been paid in the previous studies to convergence 
properties of the pairing correlation energy in finite nucleijl , to the energy dependence of the phase shift il , or to 
the renormalization with respect to the pairing gap[^|^. In the present paper we shall take a different approach 
to the cut-off, i.e., we investigate relevance of the cut-off parameter to the spatial structure of the neutron Cooper 
pair. It will be shown that the cut-off energy plays an important role to describe the strong spatial correlation at 
low density. Considering this as a physical constraint on the cut-off energy, we shall derive new parameter sets of the 
DDDI. 

Preliminary results of this work are reported in Ref. . 


II. FORMULATION 
A. BCS approximation 


We describe the neutron pair correlation in uniform neutron matter and in uniform symmetric nuclear matter 
means of the BCS approximation, which is equivalent to the Bogoliubov’s generalized mean-field approach |2; 
One of the basic equations is the gap equation, which is written in the momentum representation as 


E{k) = v'(e(fc)-M)2-t-A(fc)2. 


( 1 ) 

( 2 ) 


Here A(fc) is the pairing gap dependent on the single-particle momentum fc, while e(k) and E{k) are the single-particle 
and the quasiparticle energies. v{p — k) is the matrix element of the nucleon-nucleon interaction in the channel. 
The gap equation needs to be solved together with the number equation 


= ^ 1 
^ (27r)3 


dk 



e(fe) - 

m y’ 


(3) 


which determines the relation between the neutron density p (the Fermi momentum kp) and the chemical potential 
p. The solution of these equations defines the ground state wave function of the BCS type and the static pairing 
properties at zero temperature for a given density p. 

As the interaction acting in the channel we shall adopt a bare nucleon-nucleon force, the GSRS force[^, and the 
effective interaction given by Gogny y. The GSRS force is a local potential representation of the bare nucleon-nucleon 
interaction which is given by a superposition of three Gaussian functions: 

v{r) . (4) 


One component represents a repulsive core with the height of vi = 2000 MeV and the range parameter pi = 0.447 
fm, while two other Gaussians with V 2,3 = —240, —5 MeV and p 2,3 = 0.942, 2.5 fm represent the attraction dominant 
for 1 ^ r ^ S fm (see FigP). In spite of the simple three Gaussian representation, the GSRS reproduces rather well the 
phase shift up to about SOO MeV in the c.m. energy of the scattering nucleons. The associated scattering length 
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FIG. 1: The nucleon-nucleon potential v{r) in the channel for the GSRS and Gogny D1 forces, plotted with the solid and 
dotted curves, respectively, as a function of the relative distance r between neutrons. 


a = —17.6 fm is in close agreement with the experimental value a = —18.5 ± 0.4 fm[^. The GSRS has been used 
in some BCS calculations for the pairing at low density and for the ^P 2 pairing at high density Him. From a 
practical point of view, the analytic form makes it easy to evaluate the matrix elements of the interaction. 

The Gogny force is an effective interaction which is d esig ned for the HFB description of the pairing correlation in 
finite nuclei while keeping some aspects of the (7-matrix|52l. It is also a local potential represented as a combination 
of two Gaussian functions in the form of Eq. 0 with the range parameters /ri _2 = 0.7,1.2 fm. In the following we 
show mostly results obtained with the parameter set Dl|52| as we find no qualitative difference in the results for the 
parameter set D1S|^. A common feature of the GSRS and Gogny forces is that both are attractive in the range 
1 ^ r 3 fm while they differ largely for r 0.5 fm, where the Gogny force exhibits only a very weak repulsion instead 
of the short-range core present in the GSRS force (Fig^). The interaction range of the two forces is of the order of 
3fm. Note that the experimental effective range is Ve = 2.80 ±0.11 We shall also apply a zero-range contact 
force v{r) oc 5{r). Our treatment of this interaction will be described separately in Section IVI 

As the single-particle energy e{k) we use, in the case of the Gogny interaction, the Hartree-Fock single-particle 
spectrum derived directly from the same interaction. In the case of the bare force, it would be better from a 
viewpoint of the self-consistency to use the Brueckner Hartree-Fock spectrum. But for simplicity we adopt in the 
present analysis an effective mass approximation. Namely the single-particle energy is given by e(fc) = k'^/2m*, where 
the effective mass m* = [d^e{k)/d‘^k\kp) is derived from the Gogny HF spectrumj^ for the parameter set Dl. 

We solve the gap and number equations, 0 and 0 , without introducing any cut-off. The momentum integrations 
in the two equations are performed using a direct numerical method, where the maximum momentum kmax for the 
integration is chosen large enough so that the result does not depend on the choice of kmax- We adopt kmax = 20 fm“^ 
for the G3RS, and kmax = 10 fm“^ for the Gogny interaction. Note here that it is dangerous to introduce a small 
energy window around the chemical potential (or the Fermi energy) or a cut-off at a small momentum in evaluating 
the r.h.s. of the gap equation. Such an approximation may be justihed only in the case of the weak coupling BGS 
where the pairing gap is considerably smaller than the Fermi energy, but it is not applicable to the strong coupling 
case ^3- Note also that the chemical potential fj, and the Fermi energy ep are not the same except in the limit of the 
weak coupling. We dehne the Fermi momentum kp through the nominal relation to the density p = -^^kp, and the 
Fermi energy ep by ep = e(kp). The pairing gap /S.p = A(kp) at the Fermi momentum is used below as a measure 
of the pair correlation. In the following p denotes always the neutron density. (The total nucleon density is ptot = ‘2-p 
in the case of symmetric nuclear matter.) We define, as a reference value, the normal neutron density by po = -^kpQ 
with kpQ = 1.36 fm“^. 

In order to investigate the spatial structure of the neutron Cooper pair, it is useful to look into its wave function 
represented as a function of the relative distance between the partner neutrons of the pair. It is given by 


'i'pairir) = C {<^01 i’Hr J dkukVkP^'^'^ 


(5) 










'^k'^k 
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2E{k)’ 


( 6 ) 


in terms of the u,t!-factors except by the normalization factors C and C'. Here ^/>^(rcr) (cr =t,i) is the creation 
operator of neutron and |<i)o) is the BCS ground state. The Cooper pair wave function depends only on the relative 
distance r = \r — r'\ between the partners as it is an s-wave. We evaluate the momentum integral in Eq. © in the 
same way as in the gap and the number equations. 

It is useful to evaluate the size of the neutron Cooper pair. A straightforward measure is the r.m.s. radius of the 
Cooper pair 


where 


^rms — 


{r'^) = J drr'^\'i>pazr{r)\'^ 


ir {m^kVky 

/g°° dkk'^ {ukVkf 


(7) 

( 8 ) 


can be calculated directly from the Cooper pair wave function and/or from the u, w-factors in the momentum 

space. If one assumes weak coupling, the Pippard’s coherence length 


fi^kF 

m*TTAF 


(9) 


given analytically in terms of the gap and the Fermi momentum may be used also as another estimate of the size of 
the Cooper pair. In the following we mostly use ^rms since this quantity itself has a solid meaning even in the case 
of the strong coupling BEC case and in the crossover region between BCS and BEC. We shall use for qualitative 
discussions. 

In the present paper, we neg lect hig her order many-body effects which go beyond the BCS approximation. In many 
calculations |66l l67l 1681 l69l l70l l7lL l72l 17,1.174L I7 !t| the higher order effects in low density neutron matter are predicted 
to reduce the pairing gap by about a factor of two, which is however very much dependent on the prescriptions 
adopted 0, 0 except for the low density limit p —> o[7^. A recent Monte Carlo study [ t^ using the realistic bare 
force suggest s the gap close to the BCS result. The higher order effects in symmetric nuclear matterjT^ and in finite 
nuclei rTii 1^ are estimated to increase the gap. Keeping in mind these ambiguities, we consider that the BCS 
approximation provides a meaningful zero-th order reference. 


B. Pairing gap and coherence length 

Figure [21 shows the neutron pairing gap Ap obtained with the GSRS and Gogny D1 forces both for neutron matter 
and for symmetric matter. Results assuming the free single-particle spectrum (equivalent to the use of m* = m) 
are also shown for comparison. The pairing gap becomes maximum around p/po ~ 0.1 — 0.3 in all the cases. The 
gap decreases gradually with further decreasing the density. The difference between neutron matter and symmetric 
nuclear matter, which originates from the effective mass effect, becomes negligible at low density p/po'^5 x 10“^. 

The pairing gap obtained with the GSRS force is very similar to those obtained with more realistic models of 
the bare force (OPEGH, 113, E3|, Reid@, |4^ Argonne|4l. li^. ParisE|, Bonn^^ NijmegenQ). This is 
because the gap is essentially determined by the phase shift function]^, and the GSRS force reproduces the 
experimental phase shift, though not as accurately as the modern forces. There is small difference from them: the 
maximum gap Ap 2.5 MeV around p/po ^ 0.2 (or kF ~ 0.8fm~^) for neutron matter is slightly smaller in the 
GSRS by about 5-20%. The difference may be due to a limitation of the simple three-Gaussian representation, but 
the quality is enough for the following discussions. The gap obtained with the Gogny force is consistent with those in 
the previous calculations [d^ l5(l IhH . If we compare with the Gogny and GSRS results, they exhibit a similar overall 
density dependence, and a significant difference between the two forces is seen only at modest density p/po'^bx 10“^. 
Garrido et al. suggests that the similarity may indicate a possible cancellation among higher order effects in the 
case of symmetric matter. We find a less significant difference in the gap at rather low density p/po ^ 10“^, though 
not visible in FigEl This arises from the fact that the scattering length a = —13.5 fm of the Gogny D1 force deviates 
from the GSRS value a = —17.6 fm. 

The r.m.s. radius ^rms of the neutron Gooper pair calculated for neutron and symmetric nuclear matter with the 
bare or the Gogny forces are shown in Fig|21and Tabled The calculated ^rms is consistent with the r.m.s. radius (or 
the Pippard’s coherence length) reported in the BCS calculations using other models of the bare force pi Idil Id^ Id^ . 
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FIG. 2: The pairing gap Af in neutron and symmetric nuclear matter as a function of the neutron density p/po- The results 
for the GSRS force are shown with the symbols: cross for the free single-particle spectrum, square for neutron matter, and 
diamond for nuclear matter. The results with the Gogny D1 force are plotted with the dashed, dotted, and solid curves for 
matter with the free single-particle spectrum, for neutron matter, and for symmetric matter, respectively. 



FIG. 3: The r.m.s. radius ^rms of the neutron Cooper pair in uniform matter, plotted as a function of the neutron density 
p/po. The results for symmetric nuclear and neutron matter obtained with the Gogny D1 force are shown by the solid and 
dotted curves, respectively, while the results for symmetric nuclear and neutron matter with the GSRS force are shown by the 
diamond and square symbols, respectively. The average inter-neutron distance d = is plotted with the dot-dashed line. 


Here we would like emphasize characteristic density dependence of ^rms- It is seen from FigOthat ^rms decreases 
dramatically by nearly a factor of ten from a large value of the order of ^ 50 fm around the normal density p/po 1 
to considerably smaller values ^rms ~ 4.5 — 6 fm at density around p/po ~ 0.1. The size of the Cooper pair stays at 
small values ^rms « 5 — 6 fm in the density region p/po ~ 10“^ — 0.1. It then turns to increase, but only gradually, at 
further low density. These features are commonly seen for both neutron and symmetric nuclear matter, and for both 
the GSRS and Gogny forces. 

We would like to emphasize also the smallness of the neutron Cooper pair. This may be elucidated if we compare 
the r.m.s. radius with the average inter-neutron distance d = p“^/^ = 3.09^:^^. It is seen from FigEI that 
becomes smaller than d in a very wide range of density p/po 10“^ — 0.1. The ratio ^rms/d can reach values as 
small as ~ 0.5 at density around p/po 10“^. The relation ^rms < d, i.e. the size of the neutron Cooper pair smaller 
than the average inter-neutron distance, suggests that the neutron Cooper pair exhibits strong spatial di-neutron 
correlation. 

To understand this strong spatial correlation, it is useful to consider the ratio Ap/ep between the pairing gap 
and the Fermi energy rather than the absolute magnitude of the gap. The pairing gap Ap ft! 0.2 MeV at the 
density p/po = 1/512 for example appears small in the absolute scale, but the gap to Fermi-energy ratio amounts 
to Ap/ep ft 0.4 (see Tabled, which is larger than the value Ap/ep ft 0.25 at p/po = 1/8 where the gap is nearly 
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the maximum. If we use the Pippard’s coherence length = h'^kp/m*'KAp in place of the r.m.s. radius ^rms (this 
may be justified at least for qualitative discussion since the two quantities agree within 10-25 %, see Table 0 , the 
ratio between the r.m.s. radius and the average inter-neutron distance is related to the gap to Fermi-energy ratio 
as ^rms/d ^ ^p/d ^ 0.2eF/ Ap. Consequently we can expect in the zero-th order argument that the strong spatial 
correlation ^rmsld ^ 1 emerges when the gap to Fermi-energy ratio is larger than Ap/ep ^ 0 . 2 . This is realized in the 
present calculations in the density range p/po ^ 10 “^ — 0 . 1 . In the following we shall investigate in detail the spatial 
correlation in the neutron Cooper pair around this density range. 

We shall comment here comparison with the ^SDi neutron-proton pairing in symmetric nuclear matter. In this case 
the BCS pairing gap calculated with a realistic bare force (the Paris force) is of the order of 8 MeV at maximum [i^. 
The r.m.s. radius ^rms of the neutron-proton Cooper pair is quite small, reaching to the minimum value ^rms ~ 2 
fm at density around p/po 0.2 ^ 3 - Consequently, the r.m.s. radius ^rms becomes considerably smaller than the 
average inter-particle distance d in the density interval from p/po ^ 0.5 down to the zero density limit, where ^rms 
becomes identical to the r.m.s. radius of the deuteron^^- lo case of the neutron pairing, by contrast, the 
signature of the strong coupling ^rms ^d is obtained in the wide but limited range of the density p/po ^ 10~‘^ — 0.1. 
Apart from this difference, it is noted that the qualitative trend of the Cooper pair size, e.g. shrinking with increasing 
the density from the zero density limit, is similar to that discussed in the neutron-proton case|4(T|. 


kp 

p/po 

d 

m* /m 
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^rms 
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1.36 

1 

2.27 

0.668 

62.0 

0.64 46.60 

41.76 

0.18 

-2.91 

-2.99 

1.079 

1/2 

2.87 

0.744 

33.8 

2.03 

10.80 

9.45 

0.48 

-1.83 

-1.84 

0.68 

1/8 

4.55 

0.891 

10.9 

2.60 

4.81 

3.87 

0.81 

-0.97 

-0.90 

0.34 

1/64 

9.10 

0.984 

2.45 

0.97 

5.87 

4.71 

0.92 

-0.59 

-0.52 

0.17 

1/512 

18.20 

0.998 

0.60 

0.22 

12.05 

10.30 

0.91 

-0.62 

-0.58 

neutron matter 

, G3RS 










1.36 

1 

2.27 

0.905 

42.4 

0.14 

159.8 

144.1 

0.09 

-3.70 

-3.70 

1.079 

1/2 

2.87 

0.925 

26.1 

1.52 

11.61 

10.13 

0.47 

-1.88 

-1.85 

0.68 

1/8 

4.55 

0.969 

9.89 

2.37 

4.94 

3.90 

0.80 

-0.99 

-0.90 

0.34 

1/64 

9.10 

0.995 

2.41 

0.98 

5.90 

4.61 

0.92 

-0.60 

-0.50 

0.17 

1/512 

18.20 

0.999 

0.60 

0.24 

11.16 

9.30 

0.92 

-0.55 

-0.50 


TABLE I: The pairing gap Ap, the r.m.s. radius ^rma, the Pippard’s coherence length ^p, and the probability P{d) within 
the average inter-neutron distance d, associated with the neutron Cooper pair in symmetric nuclear matter obtained with the 
Gogny D1 interaction and in neutron matter with the GSRS force, at the neutron density p/po ~ 1,1/2,1/8, 1/64,1/512, or 
equivalently kp = 1.36,1.094,0.68,0.34,0.17 fm“^. The Fermi energy ep, the effective mass m*/m, and the average inter¬ 
neutron distance d are also shown. The parameters (l/kpa)^ and {l/kpa)A of the regularized delta interaction model are also 
listed (see the text). The units for kp and Ap,eF are fm“^ and MeV, respectively while that for d,^rma and /p is fm. 


III. SPATIAL STRUCTURE OF NEUTRON COOPER PAIR 
A. Cooper pair wave function: basics 

In examining the spatial structure of the neutron Cooper pairs, we shall focus mostly on the symmetric nuclear 
matter case obtained with the Gogny force and the neutron matter case with the GSRS force. The r.m.s. radii in 
these two cases represent a rough mean value of the four results plotted in FigEl Note also that the two r.m.s. radii 
coincide with each other within 10% in a very wide interval of density p/po = 10“^ — 0.5. It is by accident, but this 
feature can be exploited to single out influences of different interactions since the comparison can be made with the 
r.m.s. radii kept the same. 

Figures0Ka-e) show the wave function "i’pairi'f') of the neutron Cooper pair for the representative values of density 
listed in Table0 The result for neutron matter calculated with the GSRS force and that for symmetric nuclear matter 
with the Gogny D1 force are plotted in the same figure for the reasons mentioned just above. The probability density 
r^|'I'pQir,(r)p multiplied by the volume element is plotted in Fig. Etf-j). 

As a quantitative measure of the spatial correlation, we evaluate also the probability P{r) for the partners of the 
neutron Gooper pair to come close with each other within a relative distance r. It is nothing but a partial integration 
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FIG. 4: (a-e) The wave function 'i’pair{r) of the neutron Cooper pair as a function of the relative distance r between the pair 
partners at the neutron density p/po = 1,1/2,1/8, 1/64, 1/512. The solid curve is for the pair in symmetric nuclear matter 
obtained with the Gogny D1 force, while the dotted curve is for that in neutron matter with the GSRS. The vertical dotted 
line represents the r.m.s. radius ^rms of the Cooper pair in neutron matter with the GSRS while the dashed line for symmetric 
nuclear matter with Gogny Dl. Here and also in the following figures the wave function is normalized by \'i>pair (r) j^r^dr = 1. 
The arrow indicates the average inter-neutron distance d. (f-j) The same as (a-e) but for the probability density r'^\^pair 
The thin dotted line in (g) and (i) is the wave function of the fictitious ’’bound state” in the free space described in the text. 


of the probability density up to the distance r: 


/o°° \'i’pair{r^Wr^‘^dr' 


( 10 ) 


An example of this quantity is shown in Fig|3 in the case of p/po = 1/2. 

Before proceeding to the main analysis we shall first point out that the GSRS and Gogny forces provide essentially 
the same spatial structure of the Gooper pair except at very short relative distances. In Fig0I a clear difference 
between the two forces is seen at short relative distances r ^ 1 fm. (Note that the normal density case, Fig0Ka,f), is 
not relevant for this discussion since the gaps and the r.m.s. radii are very different.) Apparently the suppression of 
the wave function seen at r 1 fm in the G3RS case is caused by the strong repulsive core present in the bare force. 
The Gooper pair wave function for the Gogny force does not show this short range correlation because of the lack of 
the core. On the other hand, by looking at distances r > 1 fm slightly larger than the core radius we find that the 
Gooper pair wave functions obtained with the two forces agree quite well with each other. This observation applies 
also to the probability density r'^\'^pair{T)\'^ and the probability P(r), for which the difference at short distances 
r < 1 fm becomes barely visible as the volume element is small at such short distances. Thus the spatial structure 
of the neutron Cooper pair does not depend on whether the interaction is the bare force or the effective Gogny force. 























































FIG. 5: The probability P{r) for the partner neutrons of the Cooper pair to be correlated within a relative distance r, 
calculated at density p/po = 1/2. The result for symmetric nuclear matter with the Gogny D1 force is plotted with the solid 
curve while the dotted curve represents the result for neutron matter with the GSRS force. The dashed vertical line indicates 
the r.m.s. radius of the Cooper pair in the symmetric matter case. The vertical dotted line marks the position r = 3 fm. 


provided that two cases gives the same r.m.s. radius of the Cooper pair. In the following, we shall concentrate on 
behaviors which are common to the two interactions. 

The Cooper pair wave function '^pair{r) in the coordinate representation is reported in some of the previous BCS 
calculations adopting other models of the bare force 1^ 1^ IttL and the Gogny force ji^. Our wave function 
appears consistent with those in riefs.|4^ IS. IdM l48| although in these references the wave function is shown only 
at very limited numbers of density values and up to not very large relative distances. We found, however, that the 
shape of the Cooper pair wave function shown in R,ef.differs largely from our results (Fig^J, especially at relative 
distances smaller than several fm. 


B. Density dependence 

If we compare in Fig 01 the Cooper pair wave functions at different density values, important features show up. 
An apparent observation is that the spatial extension or the size of the Cooper pair varies strongly with the density, 
in accordance with the strong density dependence of the r.m.s radius ^rms discussed in the previous section (FiglJJ. 
We emphasize here another prominent feature. Namely the profile of the Cooper pair wave function also changes 
significantly with the density. 

At the normal density p/po = 1 (Fig01a,f)), the Cooper pair wave function is spatially extended: the r.m.s. radius 
of the Cooper pair is as large as ^rms ^ 50 fm. The profile of the Cooper pair wave function in this case exhibits 
an exponential fall-off convoluted with an oscillation. This behavior is consistent with the well known expression]^ 
r'^pairir) ^ Ko{r/Tr^p)sm{kFr) for the Cooper pair wave function in the weak coupling BCS situation. Here Kq is 
the modified Bessel function, which behaves asymptotically as Ko{r/n^p) ~ exp(—( c/tt^p)). The position 

of the first node r « approximately corresponds to the average inter-particle distance d = 3.09A:p^(= 2.3fm). 
The wave function has significant amplitude for r > d since we here have a relation ^rms,p ^ d (see Table 0. This is 
a typical behavior in the situation of the weak coupling BCS. 

The Cooper pair wave function at the density p/po = 1/8 (Fig01c,h)) is very different from that at the normal 
density. Apart from the considerably small spatial extension {^rms = 4.8 — 4.9 fm), the functional form of the wave 
function behaves quite differently. We find that amplitude of the wave function is strongly concentrated within the 
average inter-neutron distance d, and that the oscillating amplitude beyond d is quite small. This is consistent with 
the observation in the previous section that the r.m.s. radius ^rms of the Cooper pair is smaller than the average 
inter-neutron distance d in this case (cf.FigjJland Table P). The probability P{d) for the partners of the Cooper pair 
to be correlated within the inter-nucleon distance d, exceeds 0.8 (Fig|S|and Table P, indicating directly the strong 
spatial di-neutron correlation. 

The neutron Cooper pair wave functions at p/po = 1/64 and 1/512, Fig0Jd,e,i,j), exhibit a behavior similar to that 
at p/po = 1/8. Inspecting more closely, we notice that the concentration within r < d is stronger than at p/po = 1/8 
while the spatial extension itself is slightly larger {^rms « 6 — 12 fm). We observe also smaller oscillating amplitude 
in the large distance region r > d (Fig0, larger values of P(d) « 0.9, and smaller ratio ^rms/d (TableP. They all 
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point to stronger spatial di-neutron correlation at these values of density. 

The Cooper pair wave function at p/po = 1/2 (Fig01b,g)) exhibits an intermediate feature between that at the 
normal density p/po = 1 and those at p/po = 1/64— 1/512. In particular, we notice that the spatial correlation 
seen at the lower density persists to some significant extent also in this case. For example, the probability density is 
strongly concentrated to the short distance region of r^3 fm (Fig^Jg)). This is more apparent in the plot of P{r) 
shown in FigEl where we find that the probability P{r) increases steeply with increasing r from r = 0, and reaches 
~ 50% already at r = 3 fm, which is roughly the interaction range of the nucleon force. This strong concentration 
within r < 3 fm may be elucidated by comparing with what could be expected if a bound pair having the same r.m.s. 
radius {^rms = 10.8 fm in this case) was formed in the free space. (We calculate this fictitious “bound state” wave 
function by increasing the strength of the Gogny D1 potential by a numerical factor.) It is noticed that the profile of 
the Cooper pair wave function differs from the “bound state” wave function which is plotted with the thin dotted line 
in Fig0Kg). In this ’’bound state” wave function, concentration of the probability within r^3 fm is not very large, 
i.e., P(3fm) = 0.24, while the probability P(3fm) associated with the neutron Cooper pair wave function is about 
twice this value. This indicates that the spatial di-neutron correlation is also strong for the moderate low density 
region p/po 0.5. A remnant of this spatial correlation is found also at the normal density (Fig0Ka,f)), but in this 
case the concentration within the interaction range is not very large (P(3fm) = 0.21), due to the very large Cooper 
pair size ( ^rms ~ 50 fm). In contrast the Cooper pair wave function at the lower density, p/po = 1/64 (Fig^i)) for 
example, is much more similar to the ’’bound state” wave function. 

Figure El shows the overall behavior of P(3fm) and P(d) as a function of the density. It is seen that the strong 
concentration within the interaction range, say P(3fm) > 0.5, is realized in the density region p/po « 5 x 10“^ — 0.5. 
The probability P(3fm) reaches the maximum value ~ 0.7 around p/po ~ 0.1, where the r.m.s. radius is the smallest. 
At lower density p/po ^ 10“^, the probability P(3fm) decreases gradually in accordance with the gradual increase of 
the r.m.s. radius of the Cooper pair. Note however that in this density region (p/po ~ 10“'^— 10“^) the concentration 
of the probability within the average inter-neutron distance d remains very large, i.e., P{d) ^ 0.8. 

All the above analyses indicate that the spatial di-neutron correlation is strong in the quite wide density interval 
p/po ~ 10-4 - 0.5. 

It is interesting to compare our result with that in a similar analysis of the Cooper pair wave function for the ^SDi 
neutron-proton pairing. In that case the Cooper pair wave function is found to merge smoothly into the deuteron 
wave function in the low density limit j39||. Correspondingly, the r.m.s. radius of the Cooper pair approaches to that 
of the deuteron, which is much smaller than the average inter-particle dista,uce|4flj. This is interpreted as a realization 
of the BEC of the deuterons in the low density region and the BCS-BEC crossover taking place with change of the 
density m, 0,0 . In the neutron pairing case, the similarity of the Cooper pair wave function to a bound state 
wave function is found only in a limited density range p/po ~ IQ-^ — IQ-^, and it never converges to a bound state 
wave function (NB. there is no bound state in this channel in the free space). The r.m.s. radius ^rms of the Cooper 
pair is only comparable to the average inter-neutron distance d in the same density interval. Although the spatial 
correlation is strong as discussed above, these qualitative observations alone are not enough to assess whether the 
region of the BCS-BEC crossover is reached in the case of the neutron pairing. We shall investigate this issue on more 
quantitative bases in the next section. 



FIG. 6: The probability P{d) for the partner neutrons to be correlated within the average inter-neutron distance d and the 
probability P(3fm) within r = 3 fm. The solid and dotted curves are for symmetric nuclear matter obtained with the Gogny 
D1 force, while the square and diamond symbols are for neutron matter with the G3RS force. 
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IV. RELATION TO THE BCS-BEC CROSSOVER 

In the previous section, we have seen the strong spatial correlation in the neutron Cooper pair wave function at low 
density. In the present section we shall elucidate its implication by making a connection to the BCS-BEC crossover 
phenomenon. 

For this purpose, we shall first describe a reference Cooper pair wave function based on a simple solvable model of 
the BCS-BEC crossover, and then compare our results with that of the reference model. As such a reference, we adopt 
a ^5 pairing model that applies generically to a dilute gas limit of any Fermion svstems [^ 1^ . for which the 
average inter-particle distance d = is supposed to be much larger than the range of interaction. The dilute limit 
is equivalent to treating the interaction matrix elements as a constant, or assuming a contact interaction. Using the 
relation between the interaction constant and the zero-energy T-matrix or the scattering length a, the gap equation 
o is written in a regularized form: 


^ =_^ [dk(^ _LV 

47r?i^a 2(27r)3 _/ \E{k) e{k)) 


( 11 ) 


The regularized gap equation m and the number equation © are now expressed analytically in terms of some 
special functions, and are easily solvable 113 ■ In this model a dimensionless parameter 1/kpa, characterized by 
the scattering length and the Fermi momentum, is the only parameter that controls the strength of interaction, and 
hence properties of the pair correlation are determined solely by l/kpa while the length scale is given by (or the 
inter-particle distance d = 3.09fcp^), and the energy scale by the Fermi energy ep = h^kp/2m 
The gap to Fermi-energy ratio Ajep and the ratio ^rms/d between the r.m.s. radius and the average inter-particle 
distance are then monotonic functions of the interaction parameter l/fcpa l^l^ . The functional form of the Cooper 
pair wave function, defined by Eq. m, is also determined only by l/kpa, except for the length scale. Since there 
is no available analytic expression for the wave function, we evaluate Eq. © by performing the momentum integral 
numerically with a help of an explicit use of a smooth cut-off function of a Gaussian form. The cut-off scale is chosen 
large enough so that the results shown below do not depend on it. 

The range l/kpa —1 of the interaction parameter corresponds to the situation of the weak coupling BCS, for 
which the pairing gap is given by the well known formula|3 113113 A/ep ~ 8e“^ exp {n/2kpa). In the opposite range 
1/kpa ^ 1, the situation of the Bose-Einstein condensation (BEG) of bound Fermion pairs (bosons) is realized. The 
crossover between the weak coupling BCS and the strong coupling BEC corresponds to the interval —1^1/kpa^ 1, 
as described in Refs. [13113 HJj 113 ■ (The case with the infinite scattering length 1/kpa = 0 is the midway of the 
crossover, called the unitarity limit.) In the following we shall adopt 1/kpa = ±1 according to Ref.|13 as boundaries 
characterizing the crossover although the transition is smooth in nature. 

In Table im we list the values of E ms/ d and Ap/ep at the “boundaries” 1/kpa = ±1 of the crossover domain 
and at the unitarity limit 1/kpa = 0 |13ll3l- Note that the r.m.s. radius comparable to the average inter-particle 
distance 0.2 1-1, or the gap comparable to the Fermi energy 0.2 ^ Af/cf 1-3 corresponds to the BCS- 

BEC crossover domain —l'^l/kpa'^1. We have discussed in the previous section the probability P{d) for the paired 
neutrons to come closer than the average inter-neutron distance d (cf. Fig|S|l. As the same quantity is easily calculated 
also in the analytic model of the BCS-BEC crossover (the result is shown in FigjS)), this quantity may be used also as a 
measure of the crossover. The calculated boundary values corresponding to 1/kpa = ±1,0 are listed in Table UTI The 
crossover region is specified by 0.8 P{d) ^1.0 while the boundary to the strong coupling BEC regime {P{d) « 1) is 
hardly visible in this measure. 

In the case of nucleonic matter, the assumption of the dilute gas limit may be justified only at very low density 
p/Po^l0~^ (or fcF^0.05 fm“^)|3ll3|83j and hence we cannot apply the above analytic model in a direct manner 
to the region of the density p/po = 10“® — 1 which we are dealing with. In order to make the application possible, we 
shall stand on a more flexible viewpoint by regarding the interaction parameter 1/kpa as a freely adjustable variable, 
rather than by fixing it from the physical value of the neutron scattering length. We shall call the model treated in 
this way the regularized delta interaction model to distinguish from the original idea of the dilute gas limit. 

The interaction parameter 1/kpa needs to be determined then. We shall require the condition that the regularized 
delta interaction model gives, for a given value of density, the same r.m.s. radius ^rms/d as that of the microscopically 
calculated neutron Cooper pair. The parameter determined in this way may be denoted (1/kpa)^. We can also 
determine the interaction parameter to reproduce the ratio Ap/ep between the gap and the Fermi energy, which we 
shall denote (l/fcFa)A. The values of (l/fcFa){ and {l/kpa)A thus determined are listed in Tabled There is no 
sizable difference between (l/kpa)^ and {l/kpa)^- The Cooper pair wave functions obtained in this reference model 
are shown in FigEl It is hard to distinguish between the two options of 1/kpa. We now compare them with the 
neutron Cooper pair obtained with the Gogny force for symmetric matter at the three representative values of density 
p/pQ = 1,1/8 and 1/512. 
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FIG. 7: (a-c) The neutron Cooper pair wave function 'ifpairir) in the regularized delta interaction model, plotted with the 
dotted curve and the cross symbol in the cases of {l/kpo)^ and (l/A:Fffl)A, respectively. The neutron Cooper pair wave function 
in symmetric nuclear matter obtained with the Gogny D1 force is also shown by the solid curve, (d-f) The same as (a-c) but 
for the probability density pair{r)\^. 


It is seen from FigElthat the wave function of the regularized delta interaction model and and that of the neutron 
Cooper pair behave very similarly at distances far outside the interaction range, r ^ 5 fm. In contrast, we notice a 
sizable disagreement for r ^ 3 fm. The disagreement is understandable as the wave function within the interaction 
range r « 3 fm of the finite range Gogny force could not be described by the zero-range delta interaction. The 
Cooper pair wave function in the regularized delta interaction model exhibits the known divergence ’^pairi'i^) oc 1/r 
for r ^ 0 (cf. Fig[7Ka-c)), and consequently the disagreement between the Gogny model and the regularized delta 
interaction model becomes serious at very short relative distances r 1 fm. Note however that the squared wave 
function weighted with the volume element stays finite as seen FigUKd-f) and hence there is no diverging difference 
in the probability density. These observations suggest that the regularized delta interaction model can account for 
the essential features of the spatial structure of the neutron Cooper pair as far as the interaction strength 1/kpa is 
chosen appropriately. 


Ijkpa ^ 

rms/d A/ep P{d) 

-1 

1.10 

0.21 

0.807 boundary to BCS 

0 

0.36 

0.69 

0.990 unitarity limit 

1 

0.19 

1.33 

1.000 boundary to BEC 


TABLE II: The reference values of Ijkpa, ^rmsld and A/cf characterizing the BCS-BEC crossover in the regularized delta 
interaction model. 

We thus have a reference frame, i.e., the regularized delta interaction model, to which the neutron pairing is mapped. 
The question on possible relation to the BCS-BEC crossover phenomena can be addressed now quantitatively. We 
first look into the ratio ^rms/d between the r.m.s radius ^rms of the Cooper pair and the average inter-particle distance 
d. The values of ^rms/d for the neutron Cooper pair obtained with the G3RS force and the Gogny D1 interaction are 
compared in Fig^with the reference values defining the “boundaries” of the BCS-BEC crossover domain (Table ITTll. It 
is seen in Figl^that the calculated ratio ^rms/d enters the domain of the BCS-BEC crossover, 1.10 > ^rms/d{> 0.19), 
in the density interval p/po ~ 10“"^ — 0.1. Note also the calculated ratio becomes closest to the unitarity limit 
irms/d = 0.36 around the density p/po ~ 10“^. In the other way around, the weak coupling BCS regime is realized 
only at very low density p/po ^ 10“'* and around the normal density p/po ^ 0.2. 

Comparing in Fig^jthe gap to Fermi-energy ratio Ap/ep with the boundary values 0.21 < Ap/ep < 1.33, we have 
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FIG. 8: The probability P{d) for the partner particles to be correlated within the average inter-particle distance d for the 
regularized delta interaction model. 


the same observation that the density region p/po ~ 10“"'^ —0.1 corresponds to the domain of the BCS-BEC crossover. 
Comparison of the third measure P{d), performed in FigEl provides us the same information. It is seen also in the 
values of (l/kpa)^ and (l/fcFa)A listed in TableOlthat the condition of the crossover region {l/kFa)^.A > —1 is met 
in the cases of p/p = 1/8,1/64 and 1/512. 

On the basis of the above analysis, we conclude that the strong spatial correlation at short relative distances seen 
in the neutron Cooper pair in the very wide density range p/pn « 10““^ — 0.1 is the behavior associated with the 
BCS-BEC crossover. 



FIG. 9: The ratio ^rms/d between the r.m.s. radius ^rms of the neutron Gooper pair and the average inter-nucleon distance d, 
calculated with the Gogny D1 force for symmetric nuclear and neutron matter (the solid and dotted lines, respectively), and 
those with the GSRS force for symmetric nuclear and neutron matter (the diamond and square symbols), plotted as a function 
of the neutron density. The reference values characterizing the BGS-BEG crossover listed in Table HH are also shown with the 
horizontal dotted and dashed lines. 


V. DENSITY-DEPENDENT DELTA INTERACTION 
A. DDDI and the cut-off energy 

The contact force whose interaction strength is chosen as a density-dependent parameter, called often the density- 
dependent delta interaction (DDDI), has been employed as a phenomenological effective interaction descri bing the 
pairing correlation in finite nuclei, especially unstable nuclei with large neutron excess flCllll|. 1^ 1^ 1^ 1^ 1^ 1^ . 
In the previous section we found that the regularized contact interaction model, which employs also the contact force 
but with an analytic regularization, describes the essential feature of the spatial structure of the neutron Cooper 
pair wave function in the whole density range. This suggests a possibility that the phenomenological DDDI may 
also describe the spatial correlation in the neutron pairing. We would like to examine from this viewpoint in what 
conditions the DDDI can be justified. 
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FIG. 10: The ratio ^pjeF between the neutron pairing gap Af and the neutron Fermi energy ew, plotted as a function of 
the neutron density. The curves and the symbol are the same as in Figl^ The reference values characterizing the BCS-BEC 
crossover listed in Table HH are also shown with the horizontal lines. 


The density-dependent delta interaction has a form 

( 12 ) 

where Vo[p] is the interaction strength which is supposed to be dependent on the density. The force acts only in the 
channel due to the projection operator (1 —Pcr)/2. It should be noted that in the applications of the density-dependent 
delta interaction to finite nuclei, an explicit and hnite cut-off energy needs to be introduced. The cut-off energy in 
this case is regarded as an additional model parameter. 

Applying the DDDI to the neutron pairing in uniform matter, the gap equation reads 

where the pairing gap here is a momentum independent constant A. For the single-particle energy e(fc), we adopt the 
effective mass approximation. We use the effective mass derived from the Hartree-Fock spectrum of the Gogny Dl. 
The momentum integral in Eq. is performed under a sharp cut-off condition 

e{k) < + Cent (14) 

where we define the cut-off energy Ccut as relative energy from the chemical potential /i. We shall treat Ccut as a 
common constant which is applied to all density. We use the same cut-off in calculating the Cooper pair wave function. 

We remark that our definition of the cut-off is different from a similar one adopted in Refs. , where a cut-off 
energy is defined with respect to the single-particle energy e(fc) measured from the bottom of the spectrum e(fc = 0), 
e.g., by imposing e(fc) < 60 MeV|58| independent of the density. In our case, by contrast, we fix an energy window 
above the chemical potential /i to Ccut- We think that the cut-off energy Ccut defined in this way can be compared 
with the quasiparticle energy cut-off Ei < Ecut adopted often in the HFB calculations for finite nuclei {Ei is the 
quasiparticle energy of the single-particle state i) mHIHIlliiElililElEI. Note that in finite nuclei 
the density varies locally with the position coordinate while the quasiparticle energy is defined globally. This may 
imply, in the sense of the local density approximation, that a fixed density-independent cut-off quasiparticle energy 
Ecut is applied to each value of local density. Note also that both our Ccut and the cut-off quasiparticle energy Ecut 
are quantities measured from the chemical potential, which approximately coincide if Ccut and Ecut are sufficiently 
larger than the pairing gap. Thus we can compare directly Ccut and Ecut, provided that Ccut is chosen as a density- 
independent constant. At the zero-density limit with /x = 0, the cut-off energy ecut simply defines an upper bound on 
the free single-particle energy e(fc) = /2m. 


B. Constraints on Ccut 

Let us investigate how the energy cut-off influences the Cooper pair wave function. We performed several calculations 
using different values of Ccut = 5,10, 30, 50,100, 200 MeV for symmetric nuclear matter. In doing so, we choose the 
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interaction strength Vb for each value of density so that the gap A calculated with a given cut-off energy Ccut coincides 
with the gap Ap obtained with the Gogny D1 force. Note that the interaction strength Vq thus determined depends 
on both the cut-off energy and the density. 


^rms [fm] 


p/po 

G-CUt “5 

10 

30 

50 

100 

200 MeV Gogny D1 

1 

48.6 

47.3 

46.7 46.6 

46.5 

46.5 

46.6 

1/2 

16.9 

13.8 

11.2 

10.9 

10.7 

10.6 

10.8 

1/8 

14.9 

10.0 

6.0 


T8 

M 

4.8 

1/64 

10.5 

7.9 

6^ 

M 



5.9 

1/512 

13.2 

12.5 

12.0 

11.9 

11.8 

11.7 

12.1 


TABLE III: The r.m.s. radius Gms of the neutron Cooper pair in symmetric nuclear matter obtained with the density dependent 
delta interaction and different cut-off energies Cent = 5,10, 30, 50,100, 200 MeV. In the rightmost column, the r.m.s. radius for 
the Gogny D1 force is also listed. The underline means that the calculated number agrees with the reference Gogny D1 result 
within 10%. 



FIG. 11: (a-c) The neutron Cooper pair wave function ’l'pair(r) in symmetric nuclear matter calculated with the DDDI having 
different cut-off energies ecut = 5,10, 30, 50,100 MeV. The result obtained with the Gogny D1 interaction is shown also by the 
solid curve, (d-f) The same as (a-c) but for the probability density r^|ff'pair(r)p. 


Table nm shows the r.m.s. radius ^rms of the neutron Cooper pair calculated with the DDDI for different values of 
the cut-off energy. Here ^rms is calculated by using the Cooper pair wave function in the coordinate representation 
evaluated up to r = 500 fm. It is seen from Table Hill that the result apparently depends on the cut-off energy Ccut- In 
the cases of p/po = 1/64 — 1/2, the dependence of ^rms on e^t is very strong. The values calculated with the small 
cut-off energy ecut = 5,10 MeV largely deviate from those obtained with the Gogny force even though the interaction 
strength is chosen to reproduce the same reference pairing gap. We consider that the small cut-off energies Ccut = 5,10 
MeV are unacceptable since they fail to describe the small size ^rms ~ 5 fm of the neutron Cooper pair at the density 
around p/po = 10“^ — 0.1. If we require that the DDDI reproduces the r.m.s. radius of the neutron Cooper pair 
within an accuracy of 10% in the whole density region of interest, use of a large value of the cut-off energy satisfying 
Ccut ^ 50 MeV is suggested. 

To see roles of the cut-off energy in more details, we show in Figll lithe neutron Cooper pair wave functions '^pairir) 
obtained for Ccut = 5,10, 30, 50,100 MeV. The plot of 'i’pair{r) indicates clearly that the Cooper pair wave function 
depends sensitively on the cut-off energy Ccut- H we adopt the small cut-off energies Ccut = 5,10 MeV the wave 
function obtained with the DDDI fails to produce the strong spatial correlation at the short relative distances r ^ 3 
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fm, which is the characteristic feature of the neutron Cooper pair wave function common to the Gogny and GSRS 
forces. (Figini shows only the Gogny result for comparison, but we remind the reader of Fig^ where the GSRS case 
is also shown.) The plot of the probability density pair{j')\'^ at the density p/po = 1/8 indicates that even the 
wave function at larger distances is not described well if the small cut-off energies Ccut = 5,10 MeV are adopted. This 
is nothing but the difficulty mentioned above in describing the r.m.s. radius with these small cut-off energies. If we 
use a large cut-off energy, say ecut ^30— 50 MeV, the wave function at large distances converges reasonably to that 
obtained with the Gogny force. Concerning the wave function at short relative distances r 3 fm, on the other hand, 
we find no convergence with respect to the cut-off energy. The value of the wave function tl>pair(0) at zero relative 
distance r = 0 increases monotonically with increasing Ccut- (Increasing further Ccut oo, will approach to 

the one for the regularized delta interaction model shown in FigEI and the value 'l'p(jj,r(0) at r = 0 will diverge.) 

It may be possible to regard ecut as a parameter which simulates the finite range of the neutron-neutron interaction. 
It is then reasonable to require that the wave function '^pair (?') of the DDDI model with an appropriate choice of ecut 
describes that of the Gogny force at distances r ^ 3 fm (within the interaction range) as well as at larger distances. In 
the case of p/po = 1/8, for example, this requirement is approximately satisfied if we choose ecut = 30 or 50 MeV, see 
Fig tTIT b'l. At p/ pq = 1/512, a good description of the wave function is obtained with ecut = 30 MeV (Fig^Jc)), and 
similarly we find ecut ~ 70 MeV for the normal density pj pq = \ fFig llir a'1'1. If we do not include in the comparison 
the wave function at very short distances r ^ 1 fm where the repulsion due to the core influences in the case of the 
bare force, the constraint on the cut-off energy may be slightly relaxed. For example, at the density p/po = 1/512, 
the wave functions for ecut = 30 and 50 MeV differ only by about ^ 20% at distances 1 < r < 3 fm, and hence the 
cut-off energy ecut = 50 MeV may also be accepted. Within this tolerance we can choose a value around ecut ^ 50 
MeV as the cut-off energy which can be used commonly in the whole density region of interest. This value can be 
compromised with the constraint ecut ^ 50 MeV which we obtained from the condition on the r.m.s. radius of the 
neutron Gooper pair. 

It is interesting to note that cut-off quasiparticle e nerg ies around Er-tp = 50 — 70 MeV have been employed in many 
of recent HFB applications to finite nuclei ]^, O, S, IH, IHE Eli ■ These cut-off energies are consistent 

with the constraint ecut ^ 50 MeV suggested from the above analysis. Much smaller cut-off energies ^ 10 MeV 
adopted in early applications of the DDDI |54l IA5| are not appropriate from the view point of the spatial structure of 
the neutron Gooper pair wave function. In Ref. the cut-off energy of 20 MeV for the single-particle energy (40 MeV 
in the center of mass frame energy) was shown to describe reasonably the scattering wave function at zero energy. 
This cut-off energy is not very different from the cut-off energy ecut ^ 30 MeV which we find most reasonable (am ong 
the selected examples) in the lowest density case p/po = 1/512. In the the delta interaction model adopted in Ref.j^ 
the cut-off energy is examined with respect to the low-energy scattering phase shift in the channel. The cut-off 
value adopted is around 5-10 MeV in the single-particle energy (9 — 20 MeV for the center of mass frame energy), 
in disagreement with our value Ccut ~ 30 MeV. The difference seems to originate from different strategies to the 
delta interaction: the momentum dependence of the interaction matrix element is accounted for by the cut-off energy 
in Ref. |59| while it is taken into account in the present approach mostly through the density dependent interaction 
strength Vo[p]. 


C. DDDI parameters 

It is useful to parameterize the interaction strength of the DDDI in terms of a simple function of the density. The 
following form is often assumed H, 0, |3 EE EE ■ 

Vo[p\ = (15) 

Pc = 0.16 fm“^ (16) 

where ptot is the total nucleon density. The parameters vo,ri and a need to be determined. In the works by Bertsch 
and Esbensen|3 and Garrido et the parameters are determined so that the parameterized DDDI reproduces the 

pairing gap obtained with the Gogny force in symmetric nuclear matter as well as the experimental s-wave scattering 
length at zero density. We shall follow a similar line, but we add the important constraint that the spatial structure 
of the neutron Gooper pair is also reasonably reproduced. As discussed in the subsection just above, this can be 
achieved if we constrain the cut-off energy to a value around ecut ^ 50 MeV. Note also our definition of the cut-off 
energy is different from that in Refs. [E EE > ^-s mentioned in Subsection IV Al 

Our procedure is as follows. We consider symmetric nuclear matter. The cut-off energy is fixed to ecut = 50 MeV 
or 60 MeV for the reasons mentioned above. We then fix the interaction strength Vo[0] = vq at zero density to a value 
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Vo which reproduces the scattering length a in the free space, vq satisfying this condition is given byUm 

(17) 

(18) 


^ kc — 7r/2a ’ 

kc 


If we use as the scattering length a in Eq. the one associated with the Gogny force, the pairing gap of the DDDI in 
the low density limit p/po ^ 0 coincides with that of the Gogny force. However, since the scattering length a = —13.5 
fm of the Gogny D1 is slightly off the experimental value, we adopt the experimental one a = —18.5 fm for Eq. C3. 
This is equivalent to constrain the DDDI at the low density limit by the bare nucleon force. 

To determine the other parameters ry and a controlling the density-dependence of the interaction strength, we first 
calculate at several representative points of density the values of Vq with which the neutron gap of the Gogny D1 
force is reproduced. We then search the parameters rj and a so that the simple function Eq. 03 fits well to the values 
of Vq thus determined. We consider the density interval p/po ^ 10“^ — 1 {kp ^ 0.3 — 1.4 fm“^) in this fitting. The 
obtained values of the parameters (denoted DDDI-Dl) are shown in Table ITVl where we list also the parameter set 
obtained when we use Ccut = 60 MeV instead of 50 MeV. Another set of the parameters derived in the same way from 
the Gogny DIS force (DDDI-DIS) is listed. We performed the same procedure also for the G3RS force (DDDI-G3RS). 

The pairing gap obtained with these parameterizations of the DDDI are shown in FiglEl The resultant gap A 
agrees with that of the corresponding reference gap to the accuracy of about one hundred keV for the whole density 
region below p/po = 1. Although the results for Ccut = 60 MeV are not shown here, the agreement with the reference 
gaps is as good as in the Ccut = 50 MeV case. 



Vo [MeV fm 

T] a 

DDDI-Dl 



Bcut = 50 MeV 

-499.9 

0.627 0.55 

Bcut = 60 MeV 

-458.4 

0.603 0.58 

DDDI-DIS 



ecut = 50 MeV 

-499.9 

0.652 0.56 

Bcut = 60 MeV 

-458.4 

0.630 0.60 

DDDI-G3RS 



ecut = 50 MeV 

-499.9 

0.872 0.58 

ecut = 60 MeV 

-458.4 

0.845 0.59 


TABLE IV: The parameter sets of the density-dependent delta interaction with the cut-off energies Ccut = 50 and 60 MeV, 
derived from the procedure applied to the Gogny D1 and DIS, and the GSRS forces. See text for details. 

It is noticed in Table HVI that the parameter a takes a similar value a = 0.58 — 0.60 (in the case of ecut = 60 MeV) for 
all of DDDI-Dl, DDDI-DIS, and DDDI-G3RS, while the difference between the Gogny forces (DDDI-D1,-D1S) and 
the G3RS force (DDDI-G3RS) is readily recognized in the value of rj. It is seen also that the difference in the cut-off 
energy influences slightly the values of uq, p and a. If we compare our result with that of Ref. |58| . our parameter 
values p = 0.60 — 0.63 for the prefactor and a = 0.55 — 0.58 for the power imply stronger density dependence in Vo[p\ 
than that in Ref. [5^ . where the parameters are determined as p = 0 . 45,0 = 0.47 from a similar fitting to the gap 
with the Gogny D1 force. This is due to the difference in the cut-off schemes mentioned in Subsection IV Al Since the 
chemical potential p, increases with the density, the energy window measured from the chemical potential p decreases 
with increasing the density in the scheme of Ref. [5^ where the cut-off e{k) < 60MeV is adopted for all density while 
in our cut-off scheme the energy window is kept constant Ccut = 50,60 MeV independent of the density. Consequently 
the stronger density dependence in Vo[p\ is needed in our case. 

It may be interesting to compare our DDDI parameter sets with those determined phenomenologically from the 
experimental pairing gap in the ground states of finite nuclei. Such a comparison is made in FigtOl where the 
density dependent interaction strength Vo[p] is plotted. The phenomenological DDDI’s employed here are the so 
called surface and mixed types, for which the prefactor parameter p is fixed to 77 = 1 and p = 0.5, respectively. The 
power parameter is assumed as a = 1 for the mixed typeHIl , while for the surface type we choose here a = 1 or 1 / 2 . 
(Note that the power parameter in the surface type DDDI was investigated in Ref. [13, and 1 / 2 ^ 0 :^ 1 is suggested 
as a reasonable range of the parameter.) The strength of the surface type DDDI used here is vq = —521 MeV fm“^ 
for a = 1 , and vq = —781 for a = 1 / 2 , taken from Ref. [T(T| . where the value of vg is determined by a Skyrme HFB 
calculation for ^^°Sn to reproduce the gap 1.25 MeV. The equivalent energy cut-off 60 MeV adopted in the HFB 
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FIG. 12: The pairing gap in symmetric nuclear matter obtained with the DDDI parameter sets shown in Table livi with the 
cut-off energy Ccut = 50 MeV. The solid, dotted and dashed curves represent the result for the parameter sets DDDI-Dl, 
DDDI-DIS, and DDDI-G3RS, respectively. The symbols represent the gap for the reference calculations with the Gogny D1 
(square) and DIS (cross) forces, and the GSRS force (circle). 



FIG. 13: The density-dependent interaction strength Vo[p] of the DDDI for the parameter sets DDDI-Dl (solid curve) and 
DDDI-G3RS (dotted curve) with ecut = 60 MeV. For comparison, Vb[p] for the phenomenological DDDI parameters of the 
surface and mixed types are also shown by the symbols. See the text for details. 


calculation corresponds to our cut-off Ccut = 60 MeV. In the case of the mixed type DDDI, the adopted strength is 
vq = —290 MeV fm“^ derived from the same condition on ^^°Sn. It is seen in Figna that the density dependence in 
the present parameterizations of Vo[p] is mild in the range ptot/pc^ 0.4, resembling more to that of the mixed type 
DDDI than to that of the surface type. At lower density ptot!Pc ^ 0.3 the density dependence is stronger than that of 
the mixed type DDDI. Note that the interaction strength Vb[p] itself is not large: it takes the values between those of 
the mixed and the surface type DDDI’s in this low density region. The above comparisons suggest that the density 
dependence in the present parameterization of Vo[p] may not be very unrealistic for applications to finite nuclei. It 
could be also suggested that the present parameterization will be free from the problems pointed to in Ref. M for 
strongly density dependent DDDI’s such as the surface type DDDI with small values of the power a ^ 1/2. For more 
definite conclusions, we need to make more quantitative analyses using the HFB calculation performed directly for 
finite nuclei. It is also interesting to compare with a new approach to the DDDI with a microscopically derived cut-off 
factors}^. This is, however, beyond the scope of this paper, and such analyses will be pursued in future. 









18 


VI. CONCLUSIONS 

We have analyzed the spatial structure of the neutron Cooper pair obtained with the BCS approximation for 
neutron and symmetric nuclear matter using a bare force and the effective Gogny interaction. The size of the Cooper 
pair varies significantly with the density: its r.m.s. radius ^rms becomes as small as ~ 5 fm around p/po « 10“^ — 0.2, 
and ^rms smaller than the average inter-neutron distance d is realized in a very wide range p/po ~ 10“^ — 0.1 at 
low density. The analysis of the Cooper pair wave function indicates that the probability for the spin up and down 
neutrons in the pair to be correlated within the average inter-neutron distance d exceeds more than 0.8 in this density 
range. The strong spatial correlation at short relative distances is also seen for modest density p/po ~ 0.5, at which 
the concentration of the pair neutrons within the interaction range ~ 3 fm reaches about 0.5. These observations 
suggest that the spatial di-neutron correlation is strong, at least in the level of the mean-field approximation, in 
low-density superfluid uniform matter in the wide range of density p/po ~ 10“^ — 0.5. The essential feature does not 
depend on the interactions. 

We have investigated the behaviors of the strong di-neutron correlation in connection with the crossover phenomenon 
between the conventional pairing of the weak coupling BCS type and the Bose-Einstein condensation of the bound 
neutron pairs. Comparing with the analytic BCS-BEC crossover model assuming a contact interaction, we found that 
the density region p/po « 10“"^ — 10“^ corresponds to the domain of the BCS-BEC crossover. 

We have examined also how the density dependent delta interaction (DDDI) combined with a finite cut-off energy 
can describe the spatial correlation of the neutron Cooper pair. The spatial correlation at short relative distances and 
the r.m.s. radius of the pair are described consistently in a wide density region 0 < p/po ^ 1 provided that we adopt 
a cut-off energy around Ccut ~ 50 MeV dehned with respect to the chemical potential. We have derived a possible 
parametrization of the DDDI, which satishes this new condition on top of the constraints on the gap in symmetric 
nuclear matter and on the scattering length in the free space. The new DDDI parameterizations may be consistent 
with or at least not strongly contradictory to the phenomenological DDDI’s derived from the gap in finite nuclei. 
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